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Abstract 

We analyze the famous Wolf's sunspot numbers. We discovered that the distribu- 
tion of the sunspot number fluctuations is a mixture of the BHP distribution with 
the Gaussian distribution. 
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1 Introduction 



The sunspots are relatively dark areas on the surface of the sun caused by 
strong concentration of magnetic flux. Lu et al. [7] stated that the physical 
picture that arises, associated to sunspots, is that solar flares are avalanches 
of many reconnection events analogous to avalanches of sands in the models 
of Bak et al. [4]. The relation between small-scale processes and the statis- 
tics of global flare properties, that follows from the self-organized magnetic 
field configuration, provides a way to learn about the physics of unobserv- 
able small-scale reconnection processes. On the other hand, Bramwell et al. 
[2] showed a relation between self-organized criticality and the universal BHP 
distribution, named after the work of Bramwell, Holdsworth and Pinton [1]. 
We make a bridge between the two literatures, showing that the distribution 
of the sunspot numbers fluctuations is close to a mixture of the BHP with a 
Gaussian distribution. 



2 A BHP and a gaussian mixture fit of the sunspot numbers fluc- 
tuations 

The average duration of the sunspot cycle is 133 months, but cycles as short as 
9 years and as long as 14 years have been observed by Rabin et al. [8]. Following 
Bramwell et al. [3] and Janosi et al. [6], we define the sunspot numbers mean 
period w^(t) by 

1 T - 1 

w fl (t) = -Y / w(t + j*133) (1) 

1 3=0 
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and the sunspot numbers standard deviation period w a (t) by 



Wa(t) = \ W^t) , (2) 



where T = 23 is the number of observed cycles. The sunspot numbers fluctu- 
ations Wf(t) is given by 

w(t) - w^t) 

Wf{t) = w<T (t) ■ (3) 
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Fig. 1. Histogram of the sunspot numbers fluctuations with a Gaussian and BHP 
fits, in a semi- log plot. 

In Figure 1, we show a fit of the BHP pdf to the histogram of the sunspot num- 
bers fluctuations. For sunspot numbers fluctuations away of the mode value, 
the BHP pdf is close to the histogram. In the range -0.9 to -0.3, close to the 
mode value of the histogram, we notice a higher concentration of the sunspot 
numbers fluctuations compared with the BHP pdf. Hence, we add a Gaussian 
fit to describe the concentration of the sunspot numbers fluctuations in this 
range. We conclude that the histogram of the sunspot numbers fluctuations 
can be well approximated by a mixture of a Gaussian pdf with the BHP pdf. 
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